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Abstract: Large scale computer simulations are presented to investigate the 
avalanche statistics of sand piles using molecular dynamics. We could show 
that different methods of measurement lead to contradicting conclusions, pre- 
sumably due to avalanches not reaching the end of the experimental table. 

Keywords: 

The physics of an evolving sandpile has been of large interest to physicists 
and engineers and there has been done much work in this field. One of the 
most popular (or sometimes unpopular) ideas is the concept of self organized 
criticality (SOC) |l[|. It has been argued by many physicists that sandpiles can 
be described by cellular automata in two or three dimensions (e.g. [[J) and by 
stochastic cellular automata (e.g. ||) which in simulations might show SOC- 
behavior. There are many effects in nature which are supposed to reveal SOC, 
and hence there is a variety of articles investigating the theory of SOC (e.g. ||). 
When particles are dropped one after the other onto the top of a sand heap 
one observes avalanches. The time intervals in between successive avalanches 
and the size distribution of the avalanches have been of large interests to many 
scientists. Experimentalists as well as theorists investigated these quantities, 
and there is a controversy whether they obey a power law or not ||||. 

In an experimental work Jaeger et al. investigated the avalanche sizes of 
a pile contained in a box with one open side. The material flow over the 
edge of the box was measured in between a pair of capacitor plates. From 
the fluctuation of the capacity they concluded that the sizes of the avalanches 
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might not be power law distributed. They count for the size of an avalanche 
the capacity change, i.e. the mass of the particles which fall down over the edge 
of the table. Bretz et al. measured the avalanche distribution by recording the 
temporal behavior of the inclination of the heap's surface ||. Contradicting 
to 0] their results support the hypothesis of the power law distribution. 

We want to present the results of a large scale computer experiment where 
we recorded the distribution of the sizes of avalanches using both methods. 
Finally we will show that, perhaps, both measurements |7|,|8| do not contradict, 
but support each other. 



In a previous paper |9j we have shown that in simulation using two dimen- 
sional molecular dynamics of non spherical particles one can find a power law 
behavior of the size distribution of the avalanches. Our particles k have been 
built up of five spheres which are connected by springs (fig. |l|). 




Fig. 1. Each of the non spherical particles consists of five spheres. 



At rest the inner sphere touches the surrounding spheres of the same particles. 
For details of the forces acting between the particles of the same grain via 



springs and the forces acting between colliding grains see fT0|j . 



Using molecular dynamics we build up the pile by dropping the particles one 
after the other on the top of the evolving pile. A particle was dropped when 
all avalanches and fluctuations caused by the previously released grain have 
faded away, i.e. we waited until the maximum velocity of the particles comes 
very close to zero. Then the inclination of the heap was measured due to the 
following procedure: Suppose the shape of the heap of height H built up on a 
surface of width B is close to a triangle. Then its slope is given by 
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where m^ k \ x^ k ' and the mass and the position of the fcth grain and M is 



the sum of the masses of all particles M = 
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. Since our heap is close 
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to, but not an ideal triangle we calculated the height H using 

H = / h(x)dx , (2) 

where x max is the x-position of the grain which is closest to the end of the 
table. From the fluctuations of the slope due to eq. flU) we can conclude the 
approximate size of the avalanche according to a decrease in the slope: 

B 2 

A (1) M= — -p(tan$'-tan$), (3) 

where $' is the slope before the dropping event, B pa P = 30.7cm is the length 
of the table and p = 0.59 g-cm~ 2 denotes the average density of the heap. This 
method for the measurement of the size of an avalanche (indexed by A^M) 
is close to the experimental method used by Bretz et al. |§. Another method 
which was used by Jaeger et al. and by Rosendahl et al. || is to measure 
the weight of the material which reaches the end of the finite table during 
an avalanche using a balance or a capacitor. They considered the weight of 
material flowing over the border of the table to be the size of the avalanche 
A®M. Similar as in the experiment in the molecular dynamics simulation 
we calculated the mass of the particles which reach the end of the table, i.e. 
> P. 

Collecting the idle time of all computers of our department over a period of 
about one year we found enough computer power to perform a large scale 
molecular dynamics simulation. The heap was build up on a rough surface 
of width P = 30.7 cm pa 200 • ( r i /■ The radii r\ of the outer spheres of 

the particles k were equally distributed in the interval rj £ (0.05,0.11) cm 
while the radii of the inner spheres are determined by the relation r!$ = 
L( fc )/V2 - rf ) (i = 1 . . .4). is the size of the jfeth grain (s. fig. p}. In [K| 
(fig. 10) we have shown that the relation = 4 reproduces well the static 

i 

friction behavior of sandpiles. The average number of particles on the heap 
was N av = 930. 

Figure |2| shows the time series of the avalanche size due to both procedures, 
A«M and A&M. The size distributions of the data shown in fie;. |2| are drawn 
in fig. [| We find that the distribution of the avalanche sizes measured by means 
of the fluctuations in the slope A^M reveals a typical power law behavior for 
avalanches smaller than AM^ < 2 gramms, while the distribution according 
to the direct measurement of the avalanche sizes A^M does not show a power 
law behavior. We claim that the difference in both figures (0) comes from the 
fact that the second method is not able to care for those avalanches which 
do not reach the end of the table. Obviously the larger the pile the higher is 
the fraction of avalanches which do not reach the border of the table. This 
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Fig. 2. Series of 450 avalanches. Top: the avalanche size A^M concluded from the 
fluctuation in slope (equation (3)), bottom: the avalanche size A^M calculated 
from the mass of particles that reach the end of the table. The fraction of small 
avalanches is much higher for the upper figure. 
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Fig. 3. The size distribution of the avalanches in log-log-scale due to A^M (o) and 
A( 2 )M (+), respectively. The sizes are measured in gramms. The line shows the 
function freq ~ A^M^ 1 - 85 
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coincides with the observations by Jaeger et al. who found a deviation 
from the power law scaling for large systems. They found a sharply peaked 
avalanche distribution of large, system over spanning avalanches. 

For large avalanches in our simulation both types of measurements lead to 
very similar results which supports our conclusion. In agreement with the 
experimental observations by Rosendahl et al. || we find large avalanche tails 
in the distribution. In the case of large piles, the direct measurement of the 
mass fluctuations, i.e. neglecting the smaller avalanches, would lead to similar 
results as Jaeger et. al. found. 

For the waiting time distribution of the large avalanches, i.e. for the distri- 
bution of the number of dropping events in between two consecutive large 
avalanches, surprisingly we find a double peak. Fig. ^ shows the distribution 
for five different sizes of the discretization intervals. The double peak structure 
is found in all five curves, hence we assume that it is not an artifact due to the 
choice of the size of the discretization interval. So far we have no explanation 
for this behavior. 
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Fig. 4. The waiting time distribution of the intervals in between two consecutive 
large avalanches. The figure shows the number of pairs of avalanches over their time 
distance measured in dropping events. The double peak structure remains preserved 
for different discretization interval sizes. 



Although our simulation seems at least not to contradict the concept of SOC 
we should remark here that there are other serious objections against applying 
the idea of SOC in the case of sand pile avalanches (see e.g. PHI). 
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